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If Newton's constant G evolves on cosmological timescales as predicted by extended gravity theo- 
ries then Type la supernovae (SnIa) can not be treated as standard candles. The magnitude-redshift 
datasets however can still be useful. They can be used to simultaneously fit for both H{z) and G{z) 
(so that local G{z) constraints are also satisfied) in the context of appropriate parameterizations. 
Here we demonstrate how can this analysis be done by applying it to the Gold SnIa dataset. We 
compare the derived effective equation of state parameter w{z) at best fit with the corresponding 
result obtained by neglecting the evolution G{z). We show that even though the results clearly 
differ from each other, in both cases the best fit w{z) crosses the phantom divide w = —1. We then 
attempt to reconstruct a scalar tensor theory that predicts the derived best fit forms of H{z) and 
G{z). Since the best fit G{z) fixes the scalar tensor potential evolution F(z), there is no ambiguity 
in the reconstruction and the potential U{z) can be derived uniquely. The particular reconstructed 
scalar tensor theory however, involves a change of sign of the kinetic term ^'(z)'^ as in the minimally 
coupled case. 
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I. INTRODUCTION 

A diverse set of cosmological observations including the 
abundance of galaxy clusters 0, the baryon fraction in 
galaxy clusters ^J, statistics of large scale redshift sur- 
veys |3| , and the angular power spectrum of the Cosmic 
Microwave Background (CMB) [J indicate that the uni- 
verse is flat and that there is a low value of the matter 
density parameter 0.16 < f2om < 0.35. Thus, in the con- 
text of standard general relativistic cosmology there is a 
gap between fiom and D.tot = 1 required for flatness. This 
gap is usually assumed to be filled by an unknown form of 
energy called dark energy B]- In addition, the magnitude- 
redshift relation for type la supernovae (SnIa) [aiHIllj 
which can probe the recent expansion history of the uni- 
verse indicates that the universe has entered a phase of 
accelerating expansion (the scale factor obeys a > 0). 
This can be reconciled with the other cosmological data 
by assuming that the dark energy has negative pressure 
and therefore has repulsive gravitational properties (see 
for a recent review). 

The dark energy component is usually described by an 
equation of state parameter w = ^ (the ratio of the ho- 
mogeneous dark energy pressure p over the energy den- 
sity p). For cosmic acceleration, a value of w < — ^ is 
required as indicated by the Friedmann equation 
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The simplest viable example of dark energy is the cosmo- 
logical constant (w = —1). This example however even 
though consistent with present data lacks physical moti- 
vation. Questions like 'What is the origin of the cosmo- 
logical constant?' or 'Why is the cosmological constant 
10^^° times smaller than its natural scale Mpi so that 
it starts dominating at recent cosmological times (coin- 
cidence problem)?' remain unanswered. Attempts to 



replace the cosmological constant by a dynamical scalar 
field (quintessence [13) which may also couple to dark 
matter llj, have created a new problem regarding the 
initial conditions of quintessence which even though can 
be resolved in particular cases (tracker quintessence) , can 
not answer the above questions in a satisfactory way. 

An alternative approach towards understanding the 
nature of dark e nerg y is to attribute it to extensions of 
general relativity[l3 on cosmological scales. Such exten- 
sions can be expressed for example through scalar-tensor 
theories In these theories the Einstein Lagrangian of 
general relativity is replaced by a generalized Lagrangian 
of the form 



(1.2) 

where we have set 8ttG = 1 {Fq = 1) and Cm represents 
the matter fields and does not depend on $ so that the 
weak equivalence principle is satisfied. A common choice 
in the Jordan frame is to set Z — > 1 by rescaling the field 

The evolution of Newton's constant predicted in the 
context of extended gravity theories requires special care 
when comparing the predictions of these theories with 
observations This evolution induces special effects 
to the physics of SnIa [11 E IH 111. The observed 
magnitude redshift relation of SnIa can be translated to 
luminosity distance-redshift relation (which leads to the 
expansion history H{z)) only under the assumption that 
SnIa behave as standard candles. This assumption is 
justified in view of the fact that the observational light 
curves of closeby SnIa are well understood and their in- 
dividual intrinsic differences can be accounted for. Nev- 
ertheless, since the accelerated expansion of the universe 
is based on the fact that the peak luminosities of distant 
supernovae appear to be 0.20 magnitude fainter than 
predicted for an empty universe and ~ 0.25 magnitude 
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fainter than a decelerating universe with fiom = 0.3, it 
is clear that even minor unaccounted evolutionary effects 
can drastically change our current view for accelerating 
universe. The possible consequences of evolutionary ef- 
fects in Snia due to changes in the zero age mass and 
metallicity of the progenitor star have been previously 
explored Il9| who found that changes in the underlying 
population cause a change in the maximum brightness 
by about 0.1 — 0.2 magnitudes. 

The corresponding evolutionary effects due to the evo- 
lution of G in scalar tensor theories have also been stud- 
ied Ho. [l 8l. The peak luminosity of SnIa is propor- 
tional [13 to the mass of nickel synthesized which is a 
fixed fraction of the Chandrasekhar mass Mch varying 
as Mch G'^/"^. Therefore the SnIa peak luminosity 
varies like L ~ G~^l^ and the corresponding SnIa abso- 
lute magnitude evolves like 
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where the subscript denotes the local values of M and 
G. Thus, the magnitude-redshift relation of SnIa in the 
context of extended gravity theories is connected with 
the luminosity distance dj^i^z) as 

15 G(z\ 

mth{z) = Mo + UogdL{z) + -rlog^r^ (1.4) 
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In the limit of constant G this reduces to the familiar 
result. On the other hand, in scalar tensor theoriesfisj 
we have 



G(z) 1 2F -I- 4(dF/d$)2 1 
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and solar system experiments |23, |^ indicate that 



dF{z) 



— 0. Assuming flatness, the expansion 
history H{z) is obtained from 
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Therefore, by fitting rrith of (|1.4|) to the observed Gold 
SnIa dataset expressed as niobsizi) and using (|1.5|l and 
(ll.6|l we may obtain the best fit forms of both H{z) and 
G(z) assuming appropriate parameterizations. It should 
be pointed out that even though the modified magnitude- 
redshift relation H1.4|l has been known for some time |0 , 
it has not been properly utilized in studies attempting to 
constrain extended theories of gravity with SnIa data (see 
eg [23 ). This task is undertaken in what follows. 

The structure of this paper is the following: In the next 
section we use simple polynomial parameterizations of 
G(z) and H{z) to fit these functions to the Gold dataset 
using equations 11. 4|) and H1.6|l . In section III we use 
the best fits of G(z) and H{z) to construct the best fit 
scalar tensor theory ie the potential U (z) and the form 
of $'(z)^. Finally in section IV we conclude, summarize 
and discuss possible implications of our results. 



II. FITTING G{z) AND H{z) TO THE GOLD 
DATASET 



The recent expansion history of the universe is best 
probed by using a diverse set of cosmological data in- 
cluding SnIa standard candles, CMB spectrum, large 
scale structure power spectrum, weak lensing surveys etc. 
Given the present quality of cosmological data, among 
the above cosmological observations the most sensitive 
and high quality probe of the recent expansion history 
H{z) is the magnitude-redshift relation of SnIa. This 
probe will be used in the present study. 

There have been two approaches in deriving|2^ the 
functions dL{z), H{z) and w{z) from the discrete set of 
magnitude-redshift data and their errorbars. According 
to the first approach |24j. a smoothing window function 
is used to derive a continuous function m{z) (or equiv- 
alently dL{z)) assuming no evolution of G in equation 
(11.41) 1 and then H{z) and w{z) are obtained using the 
well known relations 1251 
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This approach works well in the context of general rel- 
ativity where G(z) = Go but it can not be used in the 
context of scalar-tensor theories (unless supplemented by 
additional cosmological observations) because equation 
(|1.4() alone can not be used to derive both c?l(2;) and 
G(z) from the single smoothed function m{z). In fact, 
even if a smoothed form of dL[z) was obtained it would 
be hard to disentangle H(z) from G{z) in equation (|1.6|l . 

The second approach [2^ |23, ll^ ll^ is based on as- 
signing particular parameterizations (which are consis- 
tent with other observations) to the function H{z) (and 
to G(z) if applicable) involving two to three parameters 
and fitting these parameters to the magnitude-redshift 
data. The goodness of fit corresponding to any set of pa- 
rameters fli, ...,a„ is determined by the probability dis- 
tribution of oi, ...,a„ i.e. 



x„) =A/-e-x'(A/,ai, 
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where 



X^(Af, ai, ...,a„) 
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{m'"''{zi) ~ m'^izi;M,ai, ...,an)y 
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m°''=i{zi) 



.(2.4) 

and A/" is a normalization factor. If prior information is 
known on some of the parameters oi, a„ then we can 
either fix the known parameters using the prior informa- 
tion or 'marginalize', i.e. average the probability distri- 
bution (|2.3|) around the known value of the parameters 
with an appropriate 'prior' probability distribution. 
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The parameters ai,...,an that minimize the 
pression 1)2. 4|l are the most probable parameter values 
(the 'best fit') and the corresponding Sn) = 

Xmin gives an indication of the quality of fit for the 
given parametrization: the smaller Xmin better the 
parametrization. The minimization with respect to the 
parameter M can be made trivially by expanding the 
of equation (|2.4|) with respect to M asj2^ 



X {ai, ..,a„) = A- 2MB + Af^C 
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A(ai, ..,a„) 
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Equation H2.5|l has a minimum for M = B/C at 



X'^(ai, a„) = A(ai, a„) - (2.7) 

Thus instead of minimizing x^(M, ai, a„) we can min- 
imize x'^(aij •••7 On) which is independent of M. Obvi- 
ously Xmin Xmin ' 

The errors are evaluated using the covariance matrix 
of the fitted parameters [s^ and the error on any cos- 
mological quantity, e.g. the equation of state w{z;pi), is 
given by: 



i=i j=i+i 



where pi are the cosmological parameters and Cij the 
covariance matrix |3lj |. 

We considered simple polynomial parameterizations 
for the functions H{z) and G{z) of the form 

H^{z) = H^{nom{l + zf + ai{l + z) + a2{l + z)^ + 
(l-Oo™-ai -aa)} (2.9) 
G(z) = Go (l + az^) (2.10) 



where the linear term in G{z) has been ignored due to 
experimental constraints on scalar tensor theories |2lj |. 
Using the parameterizations (|2.9|l and (|2.10() we have 
minimized the x'^ of (|2.7I) using the 157 datapoints of 
the Gold dataset and equations 1)1. 4|l . H1.6|) . We used a 
prior of Oom = 0.23 and we verified that our results 
are relatively insensitive to the prior of ilom used in the 
range 0.18 < flom < 0.32. The minimum was obtained 
at = 173.045 for ai = -12.35 ± 8.55, 03 = 5.41 ± 3.82 
and a = 0.05 ± 0.04. We compared our results with the 
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FIG. 1: The best fit form of H{z) in the scalar tensor (contin- 
(2.6) nous line), minimally coupled (dashed line) and ACDM (dot- 
ted) cases. The shaded region corresponds to the la region 
of the scalar-tensor best fit. 



corresponding minimally coupled case obtained by fix- 
ing a = before minimization (setting G(z) = Go at 
all times). The corresponding minimum was obtained at 
X^ = 174.168 for ai = -4.54 ± 2.52, as = 1.96 ± 1.09. 

The best fit functions for H{z), w{z) and F{z) = 
for both the scalar-tensor and minimally coupled cases 
are shown in Figs. 1, 2 and 3 respectively along with the 
1(7 (shaded) region of the scalar-tensor best fit. 

An interesting feature of Fig. 3 is the relatively large 
(about 1 5%) cosmological variation of Newton's con- 
stant. In ^3] it is shown that CMB and SDSS constraints 
allow Newton's constant to vary by a factor of 2 on cos- 
mological scales. The variation implied by our analysis 
is about 15% and is well within the cosmological con- 
straints. 

Also, it is clear that the best fit equation of state pa- 
rameter w{z) crosses the phantom divide in both the 
scalar-tensor and the minimally coupled case at about 
z ~ 0.2. This type of crossin g wh ich seems to be favored 
by the Gold Snia dataset lillillsl (but not [Mill 
by the more recent first year SNLS dataset @) has been 
the subject of extensive studies in the literature's^ as its 
reproduction is highly non-trivial in the context of most 
theoretical models[38|. 

In the scalar-tensor case w{z) does not have the usual 
meaning of the dark energy equation of state but it is 
merely defined in terms of H(z) as in equation (|2.2|l . 
In particular we use w(z) simply as an alternative way 
of plotting H{z). Such a way is useful for comparing 
with other analyses [25I |2^ based on dark energy. Those 
analyses also derive their best fit 'w{z) from their best 
fit H{z) but they simply interpret the derived best fit 
w{z) as an equation of state parameter. Therefore even 
though the interpretation is different, the comparison is 
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FIG. 2: The best fit form of w{z) in tlie scalar tensor (con- 
tinuous line) and minimally coupled (dashed line) cases. 

meaningful since the actual quantity that is compared is 
H{z). The best fit functions H{z) and F{z) will be used 
in the next section to complete the construction of the 
best fit scalar tensor theory from the Gold dataset. 



III. RECONSTRUCTING THE 
SCALAR-TENSOR LAGRANGIAN 

The derived best fit functions F{z) and H{z) may now 
be used as input [sp, 'i^l in the field equations obtained 
from the Lagrangian 11.211 in a cosmological setup to ob- 
tain the potential U{z) and the field kinetic term <i>'(z)^. 

Assuming a homogeneous $ and varying the action 
corresponding to (|1.2() in background of a flat FRW met- 
ric 

ds^ = -dt^ + a^{t){dr'^ + {dB"^ 4- sin^ B d(^^)) (3.1) 
we find the coupled system of eauations|l3j| 

ZFB'^ = p+^<^^ -3H ■ F + U (3.2) 
~2F-H = (p+p) + + F - H ■ F (3.3) 

where we have assumed the presence of a perfect fluid 
{p,p)- Eliminating <i>^ from H3.3|l . setting 

qiz) EE H{zf/Hl (3.4) 

and rescaling U ^ U ■ Hq while expressing in terms of 
redshift z we obtain 
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(3.6) 



where the prime ' denotes differentiation with respect 
to redshift (^) and we have assigned properties of matter 
{p = 0, f2o„ = %^) to the perfect fluid. 

Given the best fit form of F{z) and H{z) obtained in 
the previous section from observations, equations (|3.5() 
and H3.6|l may be used to reconstruct the U{z) and $'(z)2 
which predict the best fit forms of H{z) {q{z)) and F{z) 
for both the minimally coupled and the scalar-tensor 
cases. The resulting forms of U{z) and (f>'(z)2 are shown 
in Figs 4 and 5 respectively. 

An interesting feature of Fig. 5 is the change of sign 
of $'2(2;) that is predicted for both classes of theories (at 
best fit) at a redshift z ~ 0.2 (the same redshift where 
the phantom divide barrier is crossed). This creates a se- 
rious challenge for both classes of theories. This problem 
was well known for minimally coupled theories but it was 
believed that scalar-tensor theories which have the poten- 
tial to cross the phantom divide in a self consistent way 
|39l l4lj| could bypass this problem. Our result however 
indicates that if the Snia evolution is taken into account 



then scalar-tensor theories are faced with a similar prob- 
lem as minimally coupled models 38] in the context of 
the Gold SnIa dataset. 

As discussed in Ref. [s^, the Gold dataset mildly 
favours dynamical dark energy compared to ACDM while 
the SNLS dataset favours ACDM even in the context of 
dynamical parameterizations. Thus we anticipate that in 
the case of the SNLS dataset both of the best fit curves of 
Fig 2 (for the scalar-tensor and minimally coupled cases) 
would tend to be closer to the w=-l line. In fact, the 
minimally coupled case with the SNLS dataset has been 
discussed in Ref. Iss'l in the context of the polynomial 
parameterizations of eq. (2.9). 

IV. CONCLUSION 

We have constructed a method to utilize the magni- 
tude redshift SnIa data in the context of extended grav- 
ity theories. Assuming simple redshift parameterizations 
for H{z) and G{z) we have found their best fit forms 
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FIG. 3; The best fit form of F{z) = 1/G(z) in tlie scalar 
tensor (continuous line) and minimally coupled (dashed line) 
cases. Due to experimental constraints F{z) displays an ex- 
tremum at the present time. Consistency with the Snia data 
requires that this extremum be a maximum for F{z) (mini- 
mum for G{z)). 
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FIG. 5: The redshift evolution of the kinetic term $'(2)^ for 
the best fit forms of H{z) and F{z) for the scalar-tensor (con- 
tinuous line) and minimally coupled (dashed line) cases. The 
change of sign at a redshift z ~ 0.2 creates problems for both 
classes of theories. 
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FIG. 4: The redshift evolution of the potential for the best 
fit forms of H{z) and F{z) for the scalar-tensor (continuous 
line) and minimally coupled (dashed line) cases. 

and the corresponding error regions. The best fit form of 
G{z) indicates a slowly decreasing Newton's constant (in- 
creasing F{z) = Qjy^) at recent cosmological times. The 
corresponding best fit form of w{z) (defined through eq. 
(|2.2|) ) was found to cross the phantom divide w = — 1 for 
both a constant and a redshift dependent G. However, 
in the later case the best fit w{z) was found to vary more 
rapidly with redshift. 

The simultaneous knowledge of both F{z) and H{z) 
allows the unambiguous reconstruction of a scalar tensor 



theory by solving the generalized Friedman equations. 
The particular reconstructed scalar tensor theory using 
the Gold SnIa dataset and the specific parameterizations, 
turned out to suffer from a similar problem as the cor- 
responding minimally coupled theory (the kinetic term 
changes sign at a recent redshift). We have also consid- 
ered more complicated parameterizations for G(z) and 
H(z) which however had only a minor effect on our re- 
sults and did not seem to alleviate the above mentioned 
problem. 

Thus, if the best fit forms of H{z) and G{z) are verified 
by future SnIa datasets, this would indicate that neither 
minimally coupled nor extended quintessence are real- 
ized in Nature. In that case the SnIa data could possibly 
be consistent with either alternative extensions of gen- 
eral relativity (eg brane worlds '12]) or by a combination 
of phantom -|- quintessence scalars (quintom models42j) 
(see also |43,|4J,|43 for alternative approaches). An in- 
teresting extension of the present work would be to use 
the best fit forms of H{z) and G{z) obtained from SnIa 
and other cosmological data, in an attempt to reconstruct 
consistently alternative extended gravity theories. 

The Mathematica file with the numeri- 
cal analysis of the paper can b e found at 
'http : / /leandros . physics .uoi.gr/ snevol.html . 
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